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Quantum communication relies on optical implementations of channels, memories and repeaters. 
In the absence of perfect devices, a minimum requirement on real-world devices is that they preserve 
quantum correlations, meaning that they have some thoughput of a quantum mechanical nature. 
Previous work has verified throughput in optical devices while using minimal resources. We extend 
this approach to the quantitative regime. Our method is illustrated in a setting where the input 
consists of two coherent states while the output is measured by two homodyne measurement settings. 



I. INTRODUCTION 

The introduction of new quantum mechanical tech- 
nologies promises to fundamentally alter the way we 
communicate. Quantum key distribution (QKD), for in- 
stance, will allow us to communicate in an intrinsically 
secure way [1, 2] . But new quantum communication tech- 
nologies will require a new telecommunications infras- 
tructure, one which is quantum-enabled. That is, this 
network must be able to properly accommodate the quan- 
tum properties that quantum communications inherently 
rely on. Such a quantum network will contain many novel 
components, such as quantum memories [3], quantum re- 
peaters [4], or, most generally, quantum channels. These 
components must each operate in a strictly quantum way. 

Of course, no technology is perfect, and quantum tech- 
nologies offer a new set of practical challenges. However, 
as we have learned from QKD, perfectly ideal devices 
are not a necessity. By shifting our efforts into classical 
post-processing of data, we can deal with imperfections 
in quantum technologies. The question then becomes, 
how much imperfection can be tolerated before a device 
is no longer operating in a sufficiently quantum way? 
We can enforce a minimal quantum requirement on de- 
vices by insisting that they do not act as measure and 
prepare channels [5] (or, in the parlance of QKD, inter- 
cept and resend channels) , since communication through 
such channels is equivalent to classical communication. 
Indeed, this type of channel destroys any quantum cor- 
relations in bipartite states when one subsystem is sent 
through it. 

Of course, this is just the minimum requirement. It is 
also important to quantify the quantum behaviour, as is 
done in the field of entanglement measures, or in QKD 
through the secret key rate. For quantum channels, we 
can ask, how well does the channel preserve quantum cor- 
relations in bipartite systems, when only one subsystem 
passes through it? To study this question, we take a state 
with well-quantified quantum correlations, send one sub- 
system through the channel, and examine the output. 
We then compare the quantum correlations detectable 



in the output with the input correlations. In fact, as 
we shall see, we can test for these correlations in a so- 
called 'prepare and measure' picture, bypassing the need 
to use actual bipartite states. A strong quantum channel 
is one which preserves all or nearly all of the quantum 
correlations. This idea corresponds to what we shall call 
the quantum throughput. Such a measure would allow 
us to characterize the suitability of devices for quantum 
communication tasks. The goal of this work is to illus- 
trate that these ideas about device characterization via 
quantum throughput can be implemented in a meaning- 
ful way. Although we will make specific choices regarding 
device types or quantification measures, the basic idea re- 
mains quite general, and our scheme can be extended and 
adapted to other methods as well. 

Finally, if we picture a future quantum communica- 
tions network consisting of many components, it should 
be evident that any device-testing procedure should be 
as experimentally practical as possible. Ideally, we seek a 
testing scenario where a finite number of test states and a 
limited set of measurements are sufficient to understand 
the quantum throughput. The latter requirement is es- 
pecially important for optical systems, which are perhaps 
the most natural choice of carrier for quantum informa- 
tion. In these systems, full tomography is not really a 
practical option because of the dimension of the Hilbert 
space. We have previously examined quantum correla- 
tions in optical devices in a qualitative way [5]; in the 
present contribution, we will extend those results to pro- 
vide a quantitative picture of optical devices. 

The rest of this paper is organized as follows. In Sec. II 
we outline our quantitative device-testing scheme, focus- 
ing mainly on optical systems. We show how to estimate 
important parameters from homodyne measurements on 
the output, and how to use these estimates to make quan- 
titative statements about the optical device. In Sec. Ill, 
we give the results of this quantification procedure for a 
wide class of optical channels, and examine the strength 
of our method. Sec. IV summarizes the paper, while Ap- 
pendices A-C provide technical details and derivations. 



II. QUANTIFICATION PROCEDURE 
A. Device testing scenario 

The quantum device testing procedure we employ is 
the same as the one found in [5]. This protocol is based 
on the idea that a truly quantum channel should be dis- 
tinguishable from those channels where the input quan- 
tum state is temporarily converted to classical data be- 
fore a new quantum state is output, a so-called measure 
and prepare channel. Measure and prepare channels are 
also called entanglement-breaking channels, as the two 
notions are equivalent [6]. This provides a hint on how 
to quantify a channel's quantum throughput, namely by 
sending part of an entangled state through the channel 
and determining the amount of entanglement that still 
remains afterwards. 

To this end, imagine we have an entangled state of the 
form 
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where system A is a qubit and system B is an optical 
mode. We can assume, without loss of generality, that 
a e I, so that \a) and \—a) denote coherent states of 
opposite phase. This is an entangled state for all values 
a =/= 0, as can be seen by calculating the entropy of entan- 
glement. Keeping subsystem A isolated, an optical chan- 
nel can be probed using subsystem B of this state, fol- 
lowed by local projective measurements {|0) (0| , |1) (1|} 
by Alice and homodyne measurements {x,p} by Bob. 
These expectation values, along with the knowledge of 
Alice's reduced density matrix pa, can be used to deter- 
mine just how much of the initial state's entanglement is 
remaining. 

Of course, states like Eq. (1) may be difficult to cre- 
ate and therefore not suited for practical device testing. 
However, notice that Alice's reduced density matrix does 
not depend on what happens in the optical channel, nor 
on any of Bob's measurement results. Her expectation 
values can be completely determined from the initial 
state \i>) ab- Indeed, Alice's measurement results can 
be thought of as classical registers which merely record 
which mode state was sent through the device. This 
observation allows us to move from an entanglement- 
based (EB) picture to an equivalent 'prepare and mea- 
sure' (PM) scenario [7, 8] , in which Alice's measurements 
are absorbed into the initial state preparation. 

In a PM scenario, we retain full knowledge of pa, 
in particular the off-diagonal coherence term [pa] 01 = 
(a\—a). We must insert this additional information by 
hand into the set of expectation values for pab- This 
distinguishes the expectation values from data which 
would come from using just a classical mixture of test 
states Pab = \ (\ a ) ( a \ + I - a ) (~ a \)- Other than this, 
the procedure is the same as the EB scenario described 
above. Quantum correlations introduced in this way 
are referred to as 'effective entanglement' [21]. Using 



this convenient theoretical trick, the testing protocol can 
be accomplished simply by probing the channel using a 
source which prepares one of the two conditional states 
{| a) , |— a)} with equal probability. If the measured ex- 
pectation values, along with the inserted knowledge of 
Pa, are not compatible with any separable qubit-mode 
state, then there is (effective) entanglement and the chan- 
nel is certifiably quantum. Exploiting the duality be- 
tween the PM picture and the EB picture, we can quan- 
tify the quantum correlations remaining in the output 
state through a suitable entanglement measure. In turn, 
this can be compared to the entanglement of the state in 
Eq. (1) to determine the quantum throughput. 



B. Quantification scheme 

Our main goal in this work is to give an estimate 
of the amount of effective entanglement observable in 
an optical system after transmission through an optical 
channel. Our method is based on the following obser- 
vation: When the two, initially pure, conditional states 
{|a) , |— a)} pass through the channel, they are subject 
to loss and noise, and evolve in general to mixed states 
{po,pi} on the infinite-dimensional mode Hilbert space; 
however, since we work with coherent states, this change 
in purity comes only from the noise. Thus, for any loss 
value, if the noise introduced by the channel is not too 
high, then the output states po and p\ will still be nearly 
pure. In this case, most of the information about the 
state is still contained in a very small subspace of the full 
infinite dimensional Hilbert space. Estimating the 'most 
significant' subspace on the mode system can therefore 
be quite useful. This subspace should contain as much 
information as possible about both conditional states. 
Additionally, we will concentrate on the simplest non- 
trivial mode subspace, namely one of dimension 2. Writ- 
ing the conditional output states po and p\ in terms of 
their eigenvectors, in order of descending eigenvalues, we 
have 
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The most significant subspace is then the one formed 
using | Aq) and (A?) as basis vectors. 

Three parameters will be important to identify this 
subspace: Aq, A", and <^Aq | A^). We will estimate these 
parameters using homodyne detection. Specifically, if a is 
the annihilation operator for the mode at Bob's detector, 
then a balanced homodyne detection scheme allows us to 
measure the field quadratures, here defined as 



x = —={a) + a) , p= —={a ] -a) . (3) 

We will use the mean values (x) m , {p)o/i ana - the vari- 
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of the quadratures from both conditional states to esti- 
mate the three subspace parameters. Exactly how this 
is done will be shown in the next part. With these pa- 
rameters, we can build a 4 x 4 density matrix pp, which 
corresponds to the projection of the full qubit-mode state 
Pab onto the two-qubit subspace spanned by the basis 

{|o> a ,|i> a }<H|a°) b ,|a;) b }. 

The idea is now to bound the entanglement of the 
full state pab using the entanglement of the projection 
pp. For this, we need to exploit the strong monotonic- 
ity under local operations and classical communication 
(LOCC) property found in many entanglement measures. 
Specifically, if we perform a complete set of local mea- 
surements on a bipartite state a, which yields (perhaps 
using some classical communication) the state a m with 
probability p mi then the strong monotonicity property 
captures the idea that the entanglement should not in- 
crease, on average, under this process. In other words, 
for a given measure £, 



£{o) > y^p m £(cr m ). 



(6) 



As the name implies, this is a stronger condition than just 
monotonicity under LOCC alone. In our case, the mea- 
surement consists of projecting the mode system onto the 
most significant subspace or onto the orthogonal comple- 
ment. We denote the former projection by 
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and the latter by P 1 - = Iab — P ■ Then, if we choose an 
entanglement measure £ with the strong monotonicity 
property, we have 



£(pab) > p£ 
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where p = Tt(P PA bP) and q = Tt(P ± pabP ± )- 

For later practicality purposes, we would like to factor 
the probabilities through the entanglement measure, so 
that we work directly with unnormalized states. The un- 
normalized projected state is thus given by pp = PpabP- 
We must be careful to choose an entanglement measure 
which, in addition to being a strong monotone, can be 
defined for unnormalized states and which permits a pos- 
itive prefactor to be absorbed into the state. We will fo- 
cus on the Negativity [9-11] in this work, for which this 
choice is justified. 

We will not attempt to estimate the second term in 
Eq. (8) coming from the orthogonal projection; we only 
note that it is non-negative, so that we have the bound 



In practice, the projected matrix pp will not be fully 
characterized and will contain open parameters. On the 
other hand, some constraints can be imposed on pp from 
our knowledge of the initial conditional states and the 
homodyne measurement results, as well as natural pos- 
itivity (pp > 0) and trace constraints (Tr(pp) < 1 for 
unnormalized pp). As a final step, we must determine 
the minimal entanglement of pp compatible with all al- 
lowed values of these open parameters, subject to the 
known constraints. We are left with the final relation 



£{Pab) > min £{p P ) 

Pp 



(10) 



which will be used as the basis for calculating bounds 
on £(pab)- The next two subsections will cover how to 
estimate pp and how we minimize the entanglement over 
all compatible forms of pp. 



C. Estimating the projected state 

The first step in our method requires determining 
the projection pp of the full state pab onto the most- 
significant subspace. For this, we need to estimate the 
three parameters Ajj, A? and (Ag| A?) from the decompo- 
sition in Eq. (2). In Ref . [12], which considers the related 
problem of effective entanglement verification using het- 
erodyne measurements (i.e. full knowledge of the Q func- 
tion), several useful formulas for estimating these max- 
imal eigenvalues and overlaps are given. These bounds 
are later refined in [13], where they are used to derive 
secret key rates for continuous variable quantum key dis- 
tribution. Here, we use these bounds as a starting point 
to build up a good estimate of the projected state for 
our quantification scheme. We will roughly follow the 
notation of [13] in the following. 

First, since the conditional output states {po,Pi} have 
unit trace, their maximal eigenvalues can be parameter- 
ized by A° = 1 - ij (j = 0, 1), with i € [0, 1]. Then Eqs. 
(66) and (69) from [13] give directly the following bound: 
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£(pab) > £{pp) 



(9) 



This bound comes up several times, so it is denoted Uj 
(j = 0, 1) to make later equations more readable. Im- 
portantly, the bound can be calculated using only the 
measured variances of the conditional states. 

Estimating the overlap (Aq A°) is more involved. We 
need to derive bounds on its magnitude based on our 
available information. Again, we begin with bounds pro- 
vided in Refs. [12, 13]. With suitable relaxations, their 
bounds can be put into a specific form which will be more 
desirable for us later, as we would ultimately like to do 
a convex optimization. The specific details of this relax- 
ation are straightforward, and are outlined in Appendix 
A. We will need an additional parameter, k, which can 
be calculated directly using the measured first moments 



{(£) , (p),-}- Defining two coherent states with the same 
means as the conditional states, 
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the new parameter is given through the overlap of these 
coherent states, 



«:=| Wi>|. 



(13) 



With this definition in place, we can give the relaxed 
bounds 
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Having these bounds, obtained purely through homo- 
dyne measurements, we can now move on to estimating 
the elements of the projected density matrix pp. We can 
already estimate matrix elements of the form (A^| pj |A^) 
using Eq. (11), but to build pp we also require bounds 
on the supplementary elements (A°|pj |A°) for i ^ j. 
To get these, we first expand po into its eigenbasis, Eq. 
(2). Then, using the fact that | (A?|</>) | 2 € [0,1] for any 
normalized vector \<fr), we can easily derive the following 
bounds on the desired matrix element (see Appendix B 
for full details): 



<A;|po|A?)<(l-t/o)|<ASJ|A?)| 2 
<A?|po|A 1 )>(l-C/o)|(A«|A?)| 2 . 
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Analogous bounds can be given for (Aq pi Aq). 

Finally, we need to estimate some elements of the off- 
diagonal blocks of pp, or else there would be no way to 
differentiate an entangeld state from a classical mixture 
of the conditional states. To this end, we label the off- 
diagonal block of the full density matrix pab by poi, so 
that it is naturally split into the form 



PAB 
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Poi Pi. 



(19) 



where the diagonal blocks correspond to the two condi- 
tional states. In the PM picture, we hold full knowledge 
of the Alice's reduced density matrix 
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where c = (a\—a). Each element in Eq. (20) is the 
trace of the corresponding element in Eq. (19), so we 
can enforce the condition Tr(poi) = c - Using this as our 
starting point, and with an appropriate basis choice for 
system B, we can determine the following off-diagonal 
bounds which can be incorporated into pp: 



I (X°q\ poi |A°) I >\c\ - ^/CTb^/l - (1 - C/01 <Ag| A?) P, 

(21) 

I <A°| poi |A°) | >\c\ - y/u iy /l-(l-Uo)\(\ \\<L)\*. 

(22) 

Details on how to arrive at these inequalities can be found 
in Appendix C. 

We now have sufficient information to construct a use- 
ful estimate of the projected state. To summarize, we 
have the quantities k and Uj (j — 0,1), which can be 
calculated from measurements of the first moments and 
second moments, respectively. We want to determine pp, 
which is the projection of pab from Eq. (19) onto the 

subspace spanned by {\0) A > |1)aI ® {|^o)_b > |^i)b}- ^ e 
have estimated some of the overlaps of pab with these 
basis vectors in Eqs. (17-18) and (21-22). These esti- 
mates depend only on the input parameter c = (a\— a) 
and on the output state quantities Uq, Ui, and | (Aq A?) |. 
This last overlap quantity is itself bounded to a region 
defined by Eqs. (14-16), which depends only on Uq, Ui 
and k. Hence, for a fixed input overlap c and a fixed set 
of homodyne measurement results, we have a parameter 
region which forms a set of constraints on pp. This re- 
gion must be searched to find the minimal entanglement 
compatible with pp. We will now move on to address 
the question of how to find the minimal entanglement 
compatible with our constraints. 



D. Minimizing the entanglement 

As mentioned earlier, we will choose the Negativ- 
ity as the entanglement measure for demonstrating our 
method. In principle, we would like to find the mini- 
mal entanglement using the methods of semidefinite pro- 
gramming. But we must make some simplifications and 
relaxations which will allow us to do so. First, we exploit 
the fact that local unitary operations cannot change the 
quantity of entanglement. Therefore, without loss of gen- 
erality, we can assume that the overlap of the maximal 
cigenstates is real and positive (since this can be accom- 
plished by a relative change of phase on subsystem B). 



<AoK)>0. 
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As well, we can perform local phase changes on subsys- 
tem A, which allows us to also make the restriction 



<A°|poi|ASJ)>0. 



(24) 



The other 
(A?|poi|A?) = 



off-diagonal element of interest, 
z, is in general still a complex 



Im(z) 




Re(z) 



FIG. 1: (Color online) Relaxing the non-convex 
constraint Eq. (22) by inscribing a square within the 

circle (solid lines). Four new convex search regions 
(FIV) are defined by the intersection of each quadrant 
with the weaker constraint Eq (26). These correspond 

to the outer regions in the diagram. Better 

approximations of the circular region can be made using 

an inscribed octagon or other polygons, but this 

increases the number of convex regions that must be 

searched for the overall minimum (dotted lines). 



number. The main problem is that Eq. (22) is a 
non-convex constraint on z. To use this constraint in a 
semidefinite program, we have to replace it with a set of 
convex constraints. We accomplish this by denoting the 
right-hand side of Eq. (22) as 



|c|-VEWl-(l-0b)KA8|A2)| 2 =:r ( 25 ) 



and expanding our constraints to the region 
|Re(z)| + |Im(z)| >r. 



(26) 



This new constraint still non-convex, but we can search 
for the minimum entanglement independently in each of 
the four quadrants, where the constraints are convex (see 
Fig. 1), and take the minimum over these four searches. 
The final result will be a lower bound to the minimum 
entanglement in the region constrained by Eq. (22). 
We can extend this idea further, replacing the inscribed 
square from Fig. 1 with any other inscribed polygon. 
With more sides, we can better approximate the non- 
convex constraint Eq. (22), but this will also increase the 
number of convex subregions which must be searched to 
find the overall minimum. Numerical evidence indicates 
that the minimum entanglement is often, though not al- 
ways, found at a point outside the circle. We tested with 
an inscribed octagon and it was not found to alter the 
final results significantly. 

The final hurdle comes from the overlap (Aq I A?). Since 
the maximal eigenstates will in general have a non-zero 
overlap (indeed, for zero overlap, we will not find any en- 
tanglement in pp), we must construct an orthogonal basis 
in order to explicitly write down a matrix representing 



pp. Doing so introduces matrix elements that are both 
linear and quadratic in the overlap (Aq| A] 1 ). If the over- 
lap is used as a parameter in the semidefinite program- 
ming, this non-linear dependence becomes problematic. 
Fortunately, it turns out that to find the minimal en- 
tanglement we only need to consider the case where the 
overlap takes the largest allowed value, i.e. (A[j| A") = b u . 
The reason for this is that, for fixed values of A°, A?, 
and (AqJAj 1 ), there always exists a CPTP map on the 
B subsystem which preserves the maximal eigenvalues 
while making the corresponding overlap larger. Such a 
local map cannot increase the entanglement, so indeed 
the minimal entanglement will be found at (A[j| A") = b u . 
This useful result will be shown in detail elsewhere [14]. 



III. RESULTS 

In the previous section, we outlined a method for cal- 
culating the effective entanglement in optical systems. 
This began with the observation that we can get bounds 
just by looking at the most significant two-qubit subsys- 
tem. The remainder of Sec. II provided the necessary 
tools to allow us to calculate these bounds efficiently as 
a semidefinite program. Now that all the pieces are in 
place, we can turn to applying our scheme. 

To illustrate our quantification method, we use data 
corresponding to the action of the optical channel on the 
field quadratures, which we assume to be symmetric for 
both signal states and for both quadratures. These sym- 
metry assumptions are made solely to aid the graphical 
representation of our results, and our method does not 
rely on them. It is also important to note that, beyond 
the symmetry, we do not make any assumptions about 
how the channel works. In the absence of experimen- 
tal data, we merely parameterize the channel's effect on 
the first quadrature moments by a loss parameter and 
on the second moments by the excess noise. Specifically, 
if the means of the two conditional output states are de- 
noted by f3j from Eq. (12), then the loss is parameterized 
through the transmittivity 
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and the symmetric excess noise (expressed in Shot Noise 
Units) by 



V = 2[V j (x) 



= 2 [Vjffl- 
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(28) 



The input states are characterized entirely by the overlap 
parameter (a\—a). 

The quantification program was carried out using the 
Negativity [9-11], 
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(29) 



This measure has all the properties demanded by our 
quantification method, but more importantly, the trace 
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FIG. 2: (Color online) Minimal Negativity of pp consistent with initial overlap (a|— a) and channel loss and noise 
parameters. For comparison, the dashed lines show the Negativity for states of the form from Eq. (1), with a 

replaced by \/±a. For no loss and no noise, the bound exactly matches the initial entanglement. As either loss or 
noise increases, the bound lowers, until it becomes trivially equal to zero. Excess variances are expressed as a 

fraction of the vacuum variance (Shot Noise Units). 



norm || • ||<z> of a matrix can be computed efficiently as a 
scmidcfinite program [15]. We have normalized the Neg- 
ativity so that a maximally entangled two-qubit state has 
Af = 1. Our calculations were done in Matlab using the 
YALMIP interface [16] along with the solver SDPT3 [17]. 
Our main results are shown in Fig. 2, where the mini- 
mal Negativity of pp compatible with the initial overlap 
(a|— a) and excess noise V is given, for various values of 
the transmittivity T. This quantity gives a lower bound 
on the Negativity of the full state pab ■ The entanglement 
of the initial state, Eq. (1), is also shown as a function 
of the initial overlap in Fig. (2a). For Figs. (2b-2c), the 
modification a — > y/Ta is made to Eq. (1) for these com- 
parisons. The initial entanglement can be compared with 
the calculated bounds to help understand the quantum 
throughput of a device. In the limit of zero excess noise 
and zero loss, our entanglement bound is tight with the 
initial entanglement. 

Our bounds are quite high for very low noise, but they 
become lower as the measurement results get more noisy. 
At some point, a non-trivial entanglement bound can no 
longer be given, despite the fact that quantum correla- 
tions can still be proven for higher noise values (cf. [5]). 
As well, for larger loss values, the tolerance for excess 
noise is lower, and the region where non-trivial bounds 
can be given becomes smaller. The exact noise value 
where our bounds become trivial depends on the initial 
overlap and on the measured loss, but the highest tolera- 
ble excess noise is around 5% of the vacuum for T = 1.0. 
This shrinks to about 3% for a transmittivity of T = 0.5. 
Though the quantification region is small, it is within the 
limits of current experimental technology [18]. 

Some entanglement degradation should be expected as 
the noise is increased, but, as mentioned earlier, entan- 
glement can still be verified (though not previously quan- 



tified) under the same testing scenario up to much higher 
noise values than seen here [5, 19]. Thus, our bounds do 
not provide the full picture. The weakening of the bounds 
with higher noise is mainly due to the estimation proce- 
dure. Certain approximations become cruder (though 
still valid) as the noise increases. First, for higher noise, 
the conditional states become more mixed, spreading out 
into more of the infinite-dimensional mode Hilbert space. 
This leads to additional information being lost when we 
truncate down from pab to pp. Another problem stems 
from the bounds we use to estimate pp. Higher noise 
leads to weaker bounds on the maximal eigenvalues from 
Eq. (11), which weakens all other inequalities. 

To examine the effects of these two approximations, we 
briefly consider a simple channel where the test state, Eq. 
(1), is mixed at a 50 : 50 beam-splitter with a thermalized 
vacuum. The first moments reduce by a factor of -4= , and 
the increased variances of the output optical states can 
be determined from the mean photon number (n) of the 
thermal state. For (n) > 0, the conditional output states 
are displaced thermal states. The reason for studying this 
channel is that we can exactly determine the maximal 
eigenvalues A[j, A°, and the overlap (AqIA?). This allows 
us to study our approximations independently, since we 
decouple the effects of the two-qubit projection from the 
homodyne parameter estimation (in practice, of course, 
our quantification scheme must use both). In Fig. (3) 
we show the result of the quantification scheme, when 
this extra information is included. We see that the tol- 
erable excess noise is > 10% of the vacuum, more than 
three times what it would be if we had to estimate the 
eigenvalues and overlap using homodyne results (cf. Fig. 
(2c)). Also included in Fig. (3) is an entanglement verifi- 
cation curve, obtained using the methods of [5, 19]. Any 
points with lower noise than this verification curve must 




FIG. 3: (Color online) Results of quantification scheme 

for a simple test channel where we can exactly 
determine the two-qubit projection parameters Ag, A?, 

and (AglAi). The overall channel transmittivity is 

T=0.5, as in Fig. (2c). We also show the region where 

entanglement can be verified (points with lower noise 

than the black line). The two-qubit projection gives 

non-trivial entanglement bounds for roughly half the 

points where the entanglement can be verified. 



IV. CONCLUSION 

We have outlined a method for quantifying the effec- 
tive entanglement in qubit-mode systems using only ho- 
modyne measurement results and knowledge of the initial 
preparation. This quantification method works particu- 
larly well if the mode subsystem exhibits low noise. By 
combining this quantification scheme with a device test- 
ing scenario which uses two nonorthogonal test states, 
one can examine how strongly an optical device or exper- 
iment is operating in the quantum domain. Our scheme 
provides a useful tool for understanding the quantum na- 
ture of optical devices, especially the question of how well 
they preserve quantum correlations. 
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Appendix A: Derivation of overlap bounds 



come from entangled states. The two-qubit projection is 
tight to the entanglement verification curve for low over- 
laps. For higher values, the projection becomes weaker, 
only working to about half the noise value that the en- 
tanglement verification curve reaches. 

Ideally, we want to be able to calculate non-trivial val- 
ues for the entanglement wherever it be verified. This 
would give us a true quantitative complement to existing 
entanglement verification methods. One obvious exten- 
sion to our method would be to truncate the mode sub- 
space using the two largest eigenstates from each condi- 
tional state, or even more. In theory, this would strictly 
improve the estimates. However, in practice, this will 
increase the complexity of the quantification calculation, 
since some simplifying assumptions (i.e. certain overlaps 
are real) may no longer be valid. As well, the number of 
additional minimizations we have to do, as in our non- 
convex relaxation of Eq. (22), increases fourfold with 
each added dimension. Another approach might there- 
fore be necessary to overcome this problem. Neverthe- 
less, the quantification scheme outlined here is a useful 
method for characterizing the degree of quantumness of 
optical channels, especially when these channels intro- 
duce low noise. 



In this appendix, we derive the bounds from Eqs. (15- 
16) for the absolute value of the overlap of the maximal 
eigenstates, | (Aq| A^ 1 ) |. From [12, 13], we have the fol- 
lowing: 

Overlap Bounds. Let the largest eigenvalue of pj be 
parameterized by 



A°=:l-^, 



(Al) 



and let the fidelity between the conditional states and 
the coherent states \/3j) from Eq. (12) be given by 



and let 
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Since we cannot calculate £ or £ in practice, we now 
modify these bounds from the above form found in [12, 
13] to one involving only the parameters k (calculated 
from first moments) and the Uj (calculated from second 
moments). To do this, we make use only of the obvious 
inequality 



< ij < e 3 < U 3 . 



(A7) 



From this, we can easily derive the following auxiliary 
inequalities: 
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(A9) 



(A10) 



It is important to note that the second and third inequal- 
ities only hold so long as U < \. For symmetric noise, the 
value U = | corresponds to Vj (x) — Vj (p) = \/2 — | « 
0.914, almost twice the vacuum variance. This value is 
far outside the region where our method gives non-trivial 
bounds, so it is not an issue. Substituting the inequali- 
ties (A8-A10) into Eqs. (A5) and (A6), we arrive at the 
bounds given in Eqs. (15-16). 



Appendix B: Derivation of bounds on 
supplementary fidelities 

Here we aim to bound the quantities (A"| p 3 - |A°) for 
i =/= j, as found in Eqs. (17-18). An eigenbasis expansion 



of po leads to 

oo 

<a?|po|a?) = (i-£ )|(a°|a?)| 2 + 5:a^(a?|a^) 

OO 
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fe=i 
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A lower bound can be derived in a similar way: 
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The bounds for (Aq| /?i A§) follow by interchanging in- 
dices. 



Appendix C: Derivation of off-diagonal bounds 

This appendix outlines the derivation of the off- 
diagonal bounds from Eqs. (21-22). We completely 
know pa, which constrains that we must have Trpoi = 
(a\— a) = c. First, we consider the full density matrix 
Pab in the basis defined by (JO) , |1)} for system A and 
the eigenbasis of p , {|A§)}^L , for system B. We can 
still write this in the block form of Eq. (19), where we 
denote the diagonal elements of the block pi by {fofcl^Lo 
and the diagonal elements of the block p i by {dk}^ =Q 
(the diagonal elements of po are its eigenvalues) . Using 
the triangle inequality, we have 



\c\ = |Trpoi| 



do + y^ d h 



fe=l 



< |do|+X)|dfc 



fc=i 



From positivity of pab, we find 

CO 



(CI) 



(C2) 



fe=i 



and from the Cauchy-Schwarz inequality, 



|c| < |d, 




The first sum is just £o and the second is 1 
using the bounds from Appendix B, we get 

&o=(A l |p 1 |Ao') 

>{i-Ui)\{\l\\l)\\ 



(C4) 



which we can substitute above to obtain 

|C| < |do| + y/Foy/l-il-UJlWlXl)? 



(C5) 



Replacing do with (Aq poi Aq), we are led to the off- 
diagonal bound 

I <Ag| poi |a°) | > |c| - v^v/hm^WP. 

(C6) 
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By applying the same arguments using the eigenba- 
sis of pi, we can arrive at an analogous bound for 
|<A?| P01 |A°>|. 
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